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*]p' A generalization of the leading-order DGLAP evolution at small x is per- 

formed for the non-singlet structure function by resumming the leading-order 



DGLAP anomalous dimension to all orders in the QCD coupling. Explicit 
expressions are obtained for the non-singlet structure function of the deep 
inelastic scattering, taking into account both the double-logarithmic and the 
single-logarithmic contributions, including the running a s effects. It is shown 
that when these contributions are included, the asymptotic small-x behaviour 
is power-like, with an exponent of about 0.4. 
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I. INTRODUCTION 

Non-singlet structure functions, i.e. flavour- dependent contributions to the deep inelastic 
structure functions, have been the object of intensive theoretical investigation. Besides being 
interesting in themselves, they are also comparatively technically simple for analysis, and 
can be regarded as a starting ground for a theoretical description of DIS structure functions. 
The standard and the most widespread way for theoretical investigation of DIS structure 
functions is the DGLAP evolution equations |T[j. As is well known, in the description of 
DIS, these equations provide also good agreement with experimental data when they 
are used in the small- x region, although they are not expected to work very well in this 
kinematical region. Actually, in the small-x region, a Regge-type behaviour for the DIS 
structure functions cannot be obtained from DGLAP, which predicts instead a dependence 
of the type f NS = ev ln ^ 1//:I ' )ln ^ 2 ^^. In perturbative QCD, such Regge-like behaviour was 
obtained |3j a long time ago for structure functions of unpolarized DIS and recently j|,[5| 
for structure functions of the polarized DIS. However, all those results are obtained in the 
leading logarithmic approximation, which has the unsatisfactory feature that a s is regarded 
as a constant. 

In the present paper in addition to providing formulae for the asymptotics of the non- 
singlet structure function at very small x, we present also an explicit expression for fj^s for 
small but finite values of x, which accounts for both leading (double-logarithmic) and sub- 
leading (single-logarithmic) contributions to all orders in QCD coupling, including running 
a s effects. As logarithms of x and of Q 2 are important at small x, we account for them both. 
In order to take running a s effects into account, we use in the present work the approach 
of ||, with some important changes. Earlier, the running a s effects were incorporated into 
the double-logarithmic equations for the non-singlet structure function at small x in HH, 
but other single logarithmic contributions were not taken into account in those works. In 
the present work we take them into account, which makes our approach self-consistent. 

In the DGLAP equations, the expressions for the anomalous dimensions take into ac- 



count a finite number of NLO contributions; in this approach, the QCD coupling depends 
on Q 2 . After the total resummation of the singular NLO contributions to the anomalous di- 
mensions with the double logarithmic accuracy performed in |||| , there appeared works on 
the small- x behaviour of the structure functions combining, to a certain extent, the results 
of the DGLAP with the results of [fH|. In particular, the works []TU|-|TB"| treat the running 
QCD coupling at small x, allowing a s to depend on the incoming photon virtuality Q 2 . We 
arrive at a different conclusion. 

The paper is organized as follows: In Sect. 2 we derive the small-x evolution equations 
f° r Ins an d solve them. In Sect. 3 we obtain the asymptotic behaviour of /ats and discuss 
it. Finally, Sect. 4 contains our conclusions. 

II. SMALL-X EVOLUTION EQUATION FOR NON-SINGLET STRUCTURE 

FUNCTIONS 

In the Born approximation, and suppressing unimportant factors, we define the non- 
singlet structure function as 

fSr=so.-x). (i) 

For the radiative corrections to the Born approximation in the kinematical region of 
small x and large Q 2 , both logarithms of Q 2 and of x are essential and should be taken 
into account. Therefore, a two-dimensional evolution equation is needed for /jvs, combining 
the evolution with respect to both x and Q 2 . It is simpler to obtain such an equation 
for the scattering amplitude M(x, Q 2 ), with the s-channel discontinuity (imaginary part) 
proportional to Jns'- 

f NS = -% s M(x,Q 2 ), (2) 

n 

where s = (p + q) 2 = p 2 + 2pq(l — x) ~ 2pq is the Mandelstam variable corresponding to 

the squared total energy of the process under discussion. In the Born approximation, M is 



represented by the two graphs in Fig.p], although only graph (a) has a non-zero imaginary 
part in s and therefore only that graph contributes to /jvs- This is the reason why there 
is no difference in the Born approximation between the non-singlet contributions to the 
spin- dependent and to the spin-independent structure functions. Graphs contributing to 
M at higher orders in a s are obtained from both graphs in Fig.[TJ, by adding to them gluon 
propagators and quark and gluon loops. Starting from the order a 2 s1 the graphs obtained 
from Fig.|Tjb, with non-ladder gluon propagators incorporated, contribute in different ways 
to the non-singlet components of the spin- dependent and the spin-independent structure 
functions (see ||). In the present paper we consider the non-singlet contribution to the 
spin-independet structure functions F\ and F 2 . In terms of the Regge theory we consider 
the positive signature contribution to f^s- The first-loop correction to the Born graph (a) 
in Fig.[L] is represented by graphs in Fig.|2j. Having chosen the planar gauge 

d^ = g»u + — k (3) 

for the virtual gluon, we make the total contribution of all three graphs in Fig.|2] to be 
concentrated in the single graph (a). Iterating this graph, we obtain all graphs (of the 
ladder type only) contributing to fas- In order to avoid the infrared (IR) singularities 
that appear in the radiative corrections, we use Lipatov's prescription of compactifying the 
impact parameter space (see e.g. [0-U), in other words, we introduce the infrared cut-off \i 
in the transverse space (with respect to the plane formed by q and p) when integrating over 
the momenta of virtual particles. With such a cut-off acting as a mass scale, one can neglect 
quark masses and still be free from IR singularities. Until now the procedure is not different 
from the one used in the DGLAP, apart from the infrared regularization, which is a purely 
technical issue. The main difference comes when specifying the limits of integration over the 
transverse momenta of virtual particles. In the DGLAP, all ladder transverse momenta fcjj_ 
are ordered as 

h± < h ± < ... < Q\ (4) 



with the numeration running from the bottom of the ladder to the top. 

Equation (|J) shows that only the integration over ki± has fi as the lower limit. The 
lower limits for other ki± with i ^ 1 are expressed through the transverse momenta. As 
is well known, the ordering of Eq. (^) leads to considering the logarithms of Q 2 only. On 
the other hand, as we investigate the small-x region, we will lift this ordering, allowing for 
the transverse momentum of quarks in any ladder rung to reach the lowest limit k± = fi. 
Obviously, this increases the region of integration over ki±. Let us represent any graph 
contributing to M as a convolution of two amplitudes connected by the rung with minimal 
k±. In the DGLAP, such a representation is given by the graph in Fig.[| with the minimal 
k± in the lowest rung only; now, having lifted the k± ordering of Eq. (^), in the same graph, 
the minimal k± is assigned to any rung. Summing up such graphs to all loops and adding 
the Born graph, one arrives at the equation represented in Fig.f|a if we use the DGLAP 
ordering and at that of Fig.||b if one drops the k± ordering. The blobs in Fig.|] denote that 
radiative corrections to all orders in QCD coupling are taken into account. As a result of the 
infrared regularization we have chosen, the amplitude M is //-dependent now and one can 
evolute M with respect to fi. It is convenient to use the Mellin transform for the amplitude 
M: 



M( S/ y,Q7^) 



r^(^F M v^. (5) 

J-ioo Z7T2 \/i / 



As we are going to take into account both logarithms of s/fi 2 and of Q 2 / fJ 2 in every order 
in a,, we assume that 



?> 



M = M Born U(Hs/fi 2 ), ln(Q7/i 2 )), (6) 

where 

M Bom = f / ? N 

s-Q 2 + ie {) 

and U has to be properly defined. Equation (|7|) leads to Q s M Born = ti5(1 — x), which agrees 
with Eqs. (0) and (g). 



Let us note that U must have a non-zero imaginary part, Q S U, when s is positive and 
$> S U = when s is negative. This is satisfied when 



U = C/(ln(- S //0, HO 1 In) = U((Hs/^ 2 ) - iir)MQ*/f)) 



(8) 



Since ln(s//z 2 ) is assumed to be much greater than ir, one can expand Eq. (g) in the 
small-x limit as follows: 



M = (-1 + m8{l - x)) 



so that 



[/(ln(s//i i ),ln(Q7/i i ))-m 



dhi(s/ jjp) 



+ ... 



(9) 



3 S M^7T 



5(1 - x)U(\n(Q 2 /n 2 )) + 



9f/ 



<91n(s//z 2 ) 
Combining ( |I~0|) and @, we obtain, in the small-a; limit, 



(10) 



/iW " /_ 



100 ^(iXuF^QW). 



(11) 



So, once we know F we also know /jvs- Now we calculate F using the equations in Fig||. 
Differentiating the left-hand sides of the equations in Fig.^ with respect to ln/z 2 , we obtain 
that 



,dM 



dM 



f 



+ 



dM 



corresponds to 



djj? d\n(s/jj 2 ) <91n(Q 2 //i 



r, dF 
dy 



(12) 



(13) 



for the Mellin amplitude F. We have defined y = \n(Q 2 / /j, 2 ) . 

Then, differentiating with respect to ln(yU 2 ) the right-hand sides of the two equations in 
Fig.^, we arrive at the equation (the Born contribution does not depend on /i and cancels 
under the differentiation): 



I^Wy) 



57T 



;i + \lu) 



F(co,y)F B (u;) 



(14) 



if we follow the DGLAP k± ordering (HI), and at the equation 



' d_ 
dy 



+ u)F(v,y) 



1 



57T" 



;i + xu) 



F(u,y)F (uj) 



(15) 



if we do not invoke this k± ordering. 



In Eqs. ( |T1D and (|T5p we have used the notations A = 1/2 and provided the DGLAP 
amplitudes with a tilde . The expression in the squared brackets is the result of a one-loop 
integration in the right-hand side of the equations represented in Fig.|j. 

The Born DGLAP amplitude for the forward scattering of quarks 



if = 4ira s (Q 2 )C F /LU, 



(16) 



with Cf = {N 2 — 1)/2N = 4/3, whereas the amplitude Fq(u) f° r the forward scattering of 
quarks in Eq. (|TJ) is unknown yet and has to be calculated independently. 

Equation ( ]T4| ) leads to the DGLAP result for the non-singlet structure function in the 
leading order. Indeed, the solution of Eq. (|TJ]) is 



Cexp 



Q 2 dr a s (r)C F 
2 r 2tt 



00 



(17) 



and we therefore arrive at the famous leading-order DGLAP formula for f^s'- 



./. 



NS 



with A given by 



«*> duj /l\ w ~ /~rl ■ 
-— - CexpU - + A 

-joo 2m \x J \ Vuo 



j C F [Q 2 dr 



(18) 



(19) 



The quantity C is not determined a priori and has to be fixed somehow. However, if 
one assumes the delta-function input ([I]) corresponding to DIS off a quark in the Born 
approximation, one can specify C by the matching: 



M = M B 



(20) 



when y — 0, i.e. 



C = l. (21) 

Then one can finalize the leading-order DGLAP result for the non-singlet structure func- 
tion with the delta-function input ( f2"4|) as: 

1 



• /v ^Lc2^ exp(A 



+ A 



(22) 



On the other hand, in order to consider the small- a; behaviour of Jns, we have to lift the 
k± ordering of Eq. (|J) and to use Eq. (|lq) rather than Eq. (|i"4[) • The solution to (|i"5l) is 



F = Cexp ([-a; + (1 + \uo)F Q {uo)/%Tx 2 }y) , (23) 

which contains two unknown quantities, C and F , which have to be specified. 

Just as for the DGLAP, input C can be fixed by making special assumptions. Then the 



expression for f^s a ^ small x is given by the Mellin transform of Eq. (|TT|) , with the integrand 
Eq. (|23"D, provided F is known. In the case of the delta-function input, one can specify 
C. Equation (|23|) shows that C actually coincides with the Mellin amplitude M when the 
incoming photon is almost on-shell, i.e. its virtuality is /i 2 . So the equation for C is similar 
to the one for M, but it differs from ( |i"5|) on two points: first, there is no Q 2 -dependence 
in the left-hand side; second, the Born contribution to C is //-dependent. We recall that 
the Born contribution to M was independent of \i (see ([!])) and therefore it vanished under 
differentiation over //. Now we have to add the Born contribution 

1/w (24) 

to the right-hand side of the equation for C . Thus, C obeys 

C(w) = - + r^ [1/w + A]C(w)Fo(w) (25) 

with the obvious solution 

C = l/[w - (1 + Xuj)F {uj)/87t 2 ] . (26) 



So, the only ingredient that still remains unknown in the right-hand side of Eq. ( |2"3"D is 
the amplitude F Q , corresponding to the scattering of quarks. 
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The equation for Fq is shown in Fig.[5| It looks pretty similar to Eq. (|25|) for C. The 
difference is in the Born amplitude A{uj)/uj for the quark scattering, which we specify below 
in Eq. (40). Then, in analogy to Eq.([25|), we obtain 



r, / x Ma) 1 
F (u) = -^ + — 



UJ 57T 2 






foM, (27) 



which has the solution 

'u - Juj 2 - (1 + \uj)A(uj)/2it 2 



Fa = 4vr 



(28) 



1 + \u 

Finally, combining the above results, we obtain the expression for the non-singlet struc- 
ture function (cf. (|TTD): 

«^ C W ^ ex P([( 1 + A ^) F o/8vr 2 ]y) (29) 

for an arbitrary input C, or 

hs = I" £■ (i)"-#^exp (|(1 + toJFo/g^B) , (30) 

where we have used Eqs. (|26D and (p8|) and A(a;) is unspecified yet. Now let us determine 

A. In doing so we use the approach suggested in ||, with some important changes. 

Let us consider the forward quark-scattering amplitude in the Born approximation. The 
only Feynman graph with a non-zero discontinuity (imaginary part) with respect to the 
Mandelstam total energy variable s' of this process is shown in Fig.|6]. It yields 

M S = -Aiia s C F ^— , (31) 

s' + le 

with a s fixed in the Born approximation. 

Among the radiative corrections to the Born amplitude, there are those, see Fig.[7|, that 
make a s running. The contributions of the graphs of that kind lead to replacing Mq by 

B, when the QCD coupling depends on s': 

B = _47ra.(s / )C*-rT— • ( 32 ) 

s' + le 

9 



Furthermore, in the planar gauge, taking into account the radiative corrections is equiv- 
alent to iterating B (we recall that we discuss the positive signature in the present paper). 
As a result of such an iteration we arrive at the amplitude Mq(s'). In order to avoid the IR 
singularities in Mo, we again use the IR cut-off p so that s' > p? . in the denominators of 
Eqs. ( ftID and ([J2|) Equivalently, it can be expressed by replacing s' + le by s' — p 2 + le in 
the denominators in fl3l|) and (|32|). One can now apply the Mellin transform to B: 



B= J —[-) R{u). (33) 



du> ( s 

-too 2m \p 2 

In order to invert R(u) through B(s), one must use the inverse transform in Eq. (p3|). 
While the Mellin transform (|5|) is the extreme high-energy limit of the Sommerfeld - Watson 
transform, the inverse transform, which respects the analytical properties of the scattering 
amplitude under discussion, involves its discontinuities (see e.g. [0]). Then, it is appropriate 
to use this inverse transform in a form similar to the one suggested in ||. The inverse 
transform of Eq. (||) is 

1 r°° 
F(u) = -- — - dpexp(-ujp)%M(p), (34) 

sin(7ru;) Jo 

where p = ln(s/p 2 ), we will approximate sin(7ra;) ~ ttuj at small u>. Then inverting Eq. (| 
we obtain 



R(u) = - 

TIU! JO 



1 r°° 
— / dpe-" p %B(p), (35) 

KU JO 

having defined p = \n(s/p 2 ). 



Using the standard one-loop form for a s 

1 



a,(s 



bH-s/A 2 QCD ) ' 
when s > 0, we obtain 



(36) 



4C F 7T [°o 



(p + m) ( s , , 



(p + m) 2 + 7r 2 [{p + m) 2 + 7T 2 ] y/i 2 
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S(s/p 2 - 1) 



exp(-pw) , 

(37) 



where 

m = \n(fx 2 /A 2 QCD ) ; (38) 

and we assume p > Kqqd- The first term in the square brackets in Eq. fl3"7| ) corresponds to 
the imaginary part of (|36|) and the second one comes from the imaginary part of l/(s—p 2 +ie). 
We rewrite the expression for R in terms of A that we introduced in (|27| ) : 



^4 
R{u) = - (39) 



with 

m r 00 dpexp(-pu) 



At \ ACf7T 



(40) 



m? + ir 2 Jo (p + m) 2 + 7r 2 

Comparing Eq. (TO) with the inverse Mellin transform of Eq. (|3"T|), we see that A(u) 
corresponds to running a s in the Mellin space. In the next section we will discuss Eq. ( ffOp 
in some more detail. 

III. ASYMPTOTICS OF F NS 

Before making use of all results obtained in the previous section, let us dwell upon our 
expression for R and A. Equations ( P9"D and ( pUf ) differ from the similar expressions in M 
on one essential point: tt 2 was not accounted for in either || or [§,|J, though its numerical 
value is quite large with respect to other terms in the integrand and cannot be neglected a 
priori. 

As stated above, A in fl40|) stands for the QCD coupling in the Mellin space. Indeed, if 
we let a s be fixed just as in the DLA, we obtain 

Rdl = ^£i^l (41) 

instead of (|40|) . We note that in the framework of DGLAP, we would have obtained A given 
by flT9"| ) instead of A/8n 2 just as it was assumed in JTD|-4rj]. The reason of such a discrepancy 
is obvious: in the ladder Feynman graphs, the current argument for a s in the n-th ladder 

11 



rung is, actually, ft n +ikn±/ Pn 0- Only when x ~ 1 do we have (3 n+ \ ~ /3 n , and therefore all 
(3 can be dropped in the arguments of a s in every rung. Then, the DGLAP ordering (|J) 
establishes Q 2 as the upper limit of integration over k\. With the above approximations, one 
arrives in due course at flT9|). However, neither of those approximations holds in the small- x 
region, where there is no ordering in k± and the upper limit is s rather than Q 2 . Dropping 
those approximations leads to the expression for f^s at small x given by the expressions in 
Eqs. ( P9|) and ( f40|) for an arbitrary input and by Eqs. ( pPf) and (f40l) when the delta-function 



input is used. 

It is also easy to check that the corresponding DGLAP expressions can be obtained 
from Eqs. fl29|) and (|30|). Indeed, when x ~ 1, the main contribution of the Mellin factor, 
exp[u; ln(l/x)], comes from the region of rather large values of u, where u ln(l/x) < 1. That 
makes it possible to expand the exponent Fq in Eq. ( p9[) into a series in 1/u. Doing so, we 
obtain 



;i + Au;)Fo/87r 



2 



.4 



I + A 

U) 



A 



■in 2 



1 


ri i 




- + A 


UJ 


I_Cg> 



+ h-j + A +C(^ 3 ). (42) 



8vr a 

Retaining only the first term, proportional to A in the right-hand side of Eq. (f42|), and 
combining the result with Eq. (|29|) one arrives at the expression for /jvs 5 which is similar 
to the DGLAP expression in Eq. (|T8|) . In order to get a complete coincidence, yA/8ir 2 has 
to be replaced by A as given by Eq. (0). The same expansion of the integrand into a 
l/cu-series relates Eq. ( p0|) and Eq. (p2|) for /tvs with the delta-function input ([[]). 

The small- a; asymptotic behaviour of f^s can be obtained by evaluating f NS with the 
saddle-point method. As the asymptotics of the integrals (|29| ) and ( p0|) are identical, we 
consider below expression (E3), which we rewrite as follows: 



1 k n ± and /? are the standard Sudakov variables for ladder virtual quark momenta. For the Sudakov 
parametrization in the context of /ns, see e.g. @||. 
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100 duj 



Jns = C — - exp 

J -too 2m 



u>£ + |(w - ^ 2 - (1 + Xuj)A(uj)/2ti 2 ) 



(43) 



where we have used the notation £ = ln(l/x) and y = ln(Q 2 //i 2 ) • According to the saddle- 
point method, when x tends to zero, 



Ins ~ fNs(u ), 
where the stationary point ujq is the largest of the roots of 

uj£ + ^(uj - ^Juj 2 - (1 + XujA/2tT 2 )) 



(44) 



duj 



0. 



(45) 



Differentiating in (f|5|) leads to 
[uj 2 - (1 + Xuj)A/2n 2 ] ( ^ + f 



[w - AA/4tt 2 - (1 + Acj)(rfA/dcu)/(47r 2 )] 2 , (46) 



which can be solved numerically for any £, y. However, for obtaining asymptotically small- x 
behaviour a further simplification is possible: when £/y 3> 1, the right-hand side of Eq. ( J46|) 
is small with respect to the left-hand side and therefore ujq ~ uj s , where uj s is the root of 



uj 2 - (1 + Xuj)(A/2tt 2 ) = 0. 



(47) 



In other words, ujq is the position of the leading singularity (the branch point) of the 
exponent in Eq. ( f43|) in the uj plane. So, with the square root in ( f43|) being zero when 
uj = u>o, we arrive at the power-like asymptotics: 

■i2 \ w o/ 2 



Ins ~ InsM = exp [u; £ + Wo/ 2 ] 



;r 



^'0 



Q 2 

V 



(48) 



On the other hand, if one expands the exponent in fl43|) into a series in l/u> (see (|42|)) 
and accounts for an arbitrary but finite number of terms in the series, one arrives at the 
well-known DGLAP asymptotics for f^s instead of the power-like one. For example, let us 



retain only the first term in (f42|). In this case ojq is the root of 

1 



du . 



uj£ + y(A/8n 2 ) - + A 

UJ 



0. 



(49) 
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Replacing (yA/8Tr 2 ) by A{Q 2 ) (cf. (p~8|) ) in order to treat the QCD coupling as in the 
DGLAP, we obtain that the leading singularity cDq is now 



Co = VM ; (50) 

therefore, in this approximation : 

Ins ~ f N sM = e^ A ^' x \ (51) 

Thus, taking into account a finite number of NLO terms in the expressions for the 
anomalous dimensions inevitably leads to the asymptotic behaviour of the DGLAP type, 
whereas the power-like behaviour is the result of accounting for the NLO contributions to all 
orders in a s . Now let us calculate the exponent in Eq. (|48|). Using Eq. (|40|), we eventually 
obtain the equation for the position of the leading singularity uo s when x tends to zero : 

m f°° dr}t~^ 



uj 2 - (1 + Aw) ' 



0. (52) 



m 2 + it 2 Jo (rj + m) 2 + 7r 2 

Here, h is the first coefficient of the /5-function : b = (33 — 2ri/)/127r. We recall that m 
in fl5"2"D stands for ln(/i 2 / Aq CD ) and that we keep the relation fi > Aqcd throughout this 
paper. 

It is easy to solve Eq. ([52]) numerically, for any fixed set of parameters n/, /x, Aq C d- For 
example, using nf — 3, /j = 1 GeV, Aq C d = 0.1 GeV, we obtain 

cj s = 0.37. (53) 

On the other hand, the small-x behaviour of the non-singlet structure functions could 
be calculated in the double-logarithmic approximation and is given by 

Ins ~ aT WB£ (54) 

with the exponent 



cu DL = J2a? L C F /ir (55) 
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(see [|J for details). The QCD coupling af L is of course fixed in the DLA and its value 
cannot be specified within the DLA accuracy. Nevertheless, as asymptotically the DLA is 
supposed to dominate eventually over subleading contributions, we can estimate such an 
"asymptotic value" af L by identifying u s and ojdl- Doing so, we obtain 

a f L = u 2 s ti/2C f . (56) 

In particular, for the same set of parameters as given above, we obtain a® L = 0.16 which 
is the value of a s (Q0 GeV 2 ). We note that this value is not at all related to the value of 
Oi s {Q 2 ) {Q 2 being the incoming photon virtuality), which has often been used in the DLA 
calculations as an estimate for the value of a® L . We would like to further comment on the 
dependence of uj s on the approximation used in our calculation. First, if we had neglected the 
single-logarithmic contributions, putting A = and also dropping the ir 2 terms in Eq. (|52|). 
we would have obtained u s = 0.43. It corresponds to af L = 0.22 = a s (6 Gev 2 ). On the 
other hand, keeping both double- and single-logarithmic contributions (i.e. A = 1/2) and 
neglecting the 7r 2 terms, we would get u s = 0.49 and therefore, for this case, af L = 0.28 = 
a s (1.5 Gev 2 ). It is clear that the ir 2 terms, which come from analyticity, account only for a 
class of constant terms and should in principle, be neglected for consistency. However, the 
importance of keeping them into account in resummation formulae has been stressed for a 
long time ]16[. Unfortunately we have not been able to compare the higher terms in the 
expansion of /jvs Eq. (|42|) with fixed higher orders in a s calculations since the subleading 
(the single- logarithmic) contributions to /jvs, even in the first loop depend on the choice of 
the factorization and the regularization procedures. It would be very interesting to see the 
effect of SL contributions to the f^s with the negative signature (these contributions appear 
as from the order a^). At present, such calculations are still in progress. Although we have 
called Eqs. ([29]) and (|30| ) a generalization of the leading-order DGLAP, those expressions 



actually include also the singular contributions in u present in the expressions for the two- 
loop DGLAP anomalous dimension for f^s as we ll as other singular terms coming from the 
resummation of the leading-order DGLAP anomalous dimension to all orders in a s . 
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IV. CONCLUSIONS 

In conclusion, we have obtained a generalization of the leading-order DGLAP for f^g to 
the case of small x, by resumming the leading-order DGLAP anomalous dimension l/o;+l/2 
to all orders in a s and by accounting for running a s effects. In addition to the double- 
logarithmic contributions we took into account the single-logarithmic ones. Our result for 
the first-loop contribution to f^g is in agreement with [T5|], provided that identical regular- 
izations are used. We have demonstrated above that our results of Eqs. (|2"P| ) for an arbitrary 
input and (^) for the delta-function input ([!]) agree with the DGLAP when we neglect all 
terms except the first one in Eq. ([42; ). We showed that, with the single logarithmic con- 
tributions and with running a s taken into account, Jns also has the power-like asymptotic 
behaviour 

fNs-x-^{Q 2 /^r' 2 (57) 

as was obtained earlier in the DLA |4j] at asymptotically small x but with another value of 
the exponent. We have also investigated the dependence of the value of uj s on the accuracy 
of the calculations. In doing so we obtained that subleading (single- logarithmic) terms 
contribute to the asymptotic behaviour of f^s by increasing uj s in Eq. (|57j ) from its value 
u) s = 0.43 in the DLA up to uo s = 0.49. On the contrary, the 7r 2 -terms strongly suppress that 
growth, decreasing the value of u s down to u s = 0.37. This result can be regarded as an 
indication for the importance of the subleading contributions. In the present work we have 
discussed mainly the asymptotic behaviour of /ns- A detailed comparison of our results 
with the DGLAP at small but finite values of x will be given elsewhere. 
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Figure captions 

Fig.l The Born graphs for the amplitude M. 

Fig. 2 The graphs contributing to /jvs in one-loop approximations. 

Fig. 3 A typical graph for f NS in the planar gauge. 

Fig. 4 The evolution equation for M. 

Fig. 5 The evolution equation for the quark scattering amplitude. 

Fig. 6 The Born graph for the quark scattering amplitude. 

Fig. 7 The graphs contributing to the running QCD coupling. 
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FIG. 1. The Born graphs for the amplitude M. 
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FIG 2. The graphs contributing to fjys m one-loop approximations. 




FIG. 3. A typical graph for /jvs in the planar gauge. 
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FIG. 4. The evolution equation for M. 




FIG. 5. The evolution equation for the quark scattering amplitude. 
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FIG. 6. The Born graph for the quark scattering amplitude. 
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FIG. 7. The graphs contributing to the running QCD coupling. 



21 



